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I .  INTRODUCTION 


In  a  recent  publication,  the  absolute  rates  of  the  reactions  F  +  H2  and 
F  +  D2  were  studied  at  T  -  295  K.  The  results  for  F  +  H2  were  In  very  good 
agreement  with  earlier  studies  performed  with  a  variety  of  experimental  tech¬ 
niques.  The  lack  of  previous  absolute  rate  data  for  F  +  D2  was  noted;  how¬ 
ever,  the  ratio  of  the  rates  for  H2  and  D2  was  In  excellent  agreement  with 

earlier  relative  measurements.  These  earlier  measurements  have  been  exten- 
2-5 

slvely  reviewed. 

*^The  extension  of  the  Infrared  multlphoton  dlssoclatlon-lnfrared  chemi¬ 
luminescence  technique  to  elevated  temperatures  permitted  us  to  derive  a  set 
of  temperature-dependent  rate  coefficients  for  both  the  F  +  H2  and  F  +  D2 
chemical  reactions^  Experimental  constraints  limited  the  upper  temperature 
range  to  ~750  K;  however,  no  previous  absolute  data  have  been  reported  for 
T  >  397  K.  These  data  and  the  data  reported  in  a  concurrent  study  by  Wurzberg 
and  Houston^  at  T  =  190-373  K  provide  a  greatly  improved  data  set  for  compari¬ 
son  to  theory  and  for  modeling  of  the  HF  and  DF  chemical  lasers.  A  number  of 
classical  trajectory  calculations  on  seml-emplrlcal  potential  energy  sur¬ 
faces  have  yielded  Arrhenius  plots  for  the  F  +  H2(D2)  reactions.  Recently, 
Cohen^®  performed  extensive  parameter  variations  In  a  transition  state  theory 
description  of  the  temperature  dependence  for  these  two  processes.  These 
theoretical  treatments  provide  valuable  modeling  data,  since  real  HF  and  OF 
chemical  lasers  operate  at  elevated  temperatures  and  any  kinetic  model  that 
describes  them  must  Incorporate  accurate  T-dependent  parameters  for  the  F  -t  H2 
and  F  +  D2  reactions.  Calculations  on  electron-beam-lnltlated  H2“F2-02 


SACK  nucD 
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nlxtures  Indicia  that  move  than  60Z  of  tha'  laaar  energy  can  be  ealtted  at 
active  nadluB  tenperaturea  betwaan  500  and  2000  In  principle,  the  nodel- 

Ing  of  these  lasers  requires  not  only  thermal  rate  coefficients,  but  cross 
sections  for  translatlonally-hot  F  (and  H)  atoms, since  these  laser  systems 
arc  generally  not  buffered  by  an  Inert  gas.  Although  the  data  presented  here 
are  for  thermallzed  reaction  rates  as  a  function  of  temperature,  they  may 
guide  the  calculation  of  reaction  cross  sections  vs  relative  velocity  for 


these  two  reactions 


IT.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


The  details  of  the  apparatus  are  given  In  Ref.  1.  The  F  atoms  were 
produced  by  the  Infrared  multlphoton  dissociation  of  SF^  and  the  rate  of  the 
F  +  H2(D2)  reactions  monitored  by  the  time-dependent  chemiluminescence  from 
HF^(DF^).  Partial  pressures  of  Ar  buffer  gas,  SFg  and  H2(D2)  were  deduced 
from  molar  flow  rates  of  the  various  gases  and  the  total  system  pressure.  The 
variable-temperature  flow  cell  required  for  the  present  experiment  was  con¬ 
structed  of  45-mm  l.d.  quartz  with  Infrasil  windows  (37.5-mm  dlam)  fused  onto 
2-cm-long  side  ports  (Fig.  1).  The  cell  was  enclosed  In  a  cylindrical, 
ceramlc-f Iber „  molded  oven  (Watlow  Co.)  with  the  cell  windows  emerging  from 
the  oven  walls.  In  principle,  the  maximum  operating  temperature  for  this 
arrangement  is  ~1200  K.  The  cell  temperature  was  measured  with  a  digital 
thermocouple  (Omega  Engineering),  which  used  a  chromel-alumel  junction  encased 
In  a  stainless-steel  sheath.  At  the  Ar  buffer  pressure  of  5  Torr,  when  the 
thermocouple  probe  was  moved  across  the  diameter  of  the  flow  tube,  negligibly 
small  temperature  variations  (±3  K)  were  produced  over  the  entire  temperature 

range  studied.  The  extremely  small  gas  flows  used  prevent  serious  systematic 

13 

errors  In  temperature  determination. 

Infrared  chemiluminescence  from  the  F  +  H2(D2)  reactions  was  collected  by 
CaF2  lenses  (f/1)  arranged  as  a  compound  microscope.  After  the  radiation 
passed  through  an  appropriate  interference  filter,  It  was  refocused  onto  a 
77  K  InSb  detector  (Texas  Instruments).  The  detected  signal  was  further 
Increased  (><2)  by  the  use  of  a  spherical  A1  backing  mirror. 
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Flg»  !•  Multiphoton  dissociation— —Infrared  fluorescence  apparatus. 


The  two  Irises  Indicated  in  Fig.  1  helped  reduce  the  amount  of  blackbody 
radiation  collected  from  the  oven  assembly.  Although  the  cell  windows  re¬ 
mained  much  cooler  than  the  cell  walls,  blackbody  emission  from  the  Infrasll 
windows  prevented  measurements  at  T  >  800  K.  This  problem  was  somewhat  worse 
for  DF,  because  the  bandpass  filter  sampled  longer  wavelengths  than  the 
equivalent  filter  for  HF.  The  signal  was  processed  with  a  Blomatlon  805  tran¬ 
sient  recorder  and  accumulated  in  a  Nicolet  1072  digital  signal  averager.  A 
typical  trace  is  illustrated  in  Fig.  2.  All  data  represent  the  sum  of  eight 
laser  pulses. 

Considerable  effort  was  expended  to  explain  and  remove  the  anomalously 
large  quenching  of  HF^  and  DF^  observed  during  the  room-temperature 
studies.^  As  predicted,  a  i&ajor  fraction  of  this  quenching  resulted  from  an 
H2O  impurity  in  the  bulk  gases  and  was  effectively  removed  by  passing  the  Ar 
and  the  Ar  +  H2(D2)  mixtures  through  a  low-pressure,  glass-wool  packed  trap 
held  at  77  K. 
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t  (arbitrary  units! 


1 


I  ! 

I 

[  Fig.  2.  Time-resolved  Infrared  chemiluminescence  trace.  Channels  1 

I  I  through  512,  0.5  us/channel;  channels  513  through  1024, 

r  I  10  Ms/channel.  Partial  pressures:  SFg  *  33  mTorr,  H2  ■■ 

i  j  24.3  mTorr,  Ar  ■  3.87  Torr.  T  ■  295*K. 


J 
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III.  DATA  ANALYSIS 


In  our  earlier  study, ^  we  outlined  the  procedure  for  relating  chemi¬ 
luminescence  Intensity  from  product  HF^(DF^)  to  the  disappearance  of  the 
F-atom  reactant.  The  initial  vibrational  populations  are  produced  by 
Reactions  (1)  and  (2): 

ki (v)  ^  ^ 

F  +  H2  —  ■->  HF(v)  +  H  V  -  1,2,3  (1) 

F  +  02—^^^  DF(v)  +  D  V  *  1,...,4  (2) 

Clearly,  the  intensity  Is  given  by 

I(t)  *  Z  A^(HF(v)]  (3) 

where  A^  is  the  radiative  lifetime  of  the  vth  vibrational  level.  The  product 
vibrational  levels  are  removed  by  diffusion,  radiation,  and  impurity 
quenching,  but  primarily  these  levels  are  removed  by  the  H2  (or  D2)  that  is  In 
relatively  large  concentration: 


HF(v)  +  H2(0) 


HF(v  -  1)  +  H2(1) 


(4) 


k5(v) 


HF(v  -  1)  +  H2(0) 


(5) 


Ordinarily,  the  kinetic  system  represented  by  Eqs.  (4)  and  (5)  produces  a 
double-exponential  decay^^  of  the  fluorescence  from  HF.  Unlike  the  conven¬ 
tional  laser-induced  fluorescence  experiment  where  HF(0)  and  H2(0)  are 
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uniformly  distributed  in  the  cell,  the  present  experiment  creates  HF’ 

Initially  in  the  2  to  3-mm-dlam  focal  caustic  of  the  laser.  The  Hr(v)  and 
H2(1)  encounter  a  large  excess  of  H2(0)  as  they  rapidy  diffuse  out  of  the 
focal  volume.  Since  [H2(1)]/[H2(0)1  remains  small,  the  equilibriias  defined  by 
Eq.  (A)  is  never  established.  Therefore,  the  effective  rate  coefficient  for 
removal  of  HF(v)  by  H2(0)  in  this  experiment  is  k^(v)  +  kjCv). 

An  approximate  analytic  solution  of  the  time  derivative  of  I(t)  was 


developed  in  Ref.  1: 


kj[H2]  -  (k^(l)(H2]  +  C2) 


{exp  [-(k.(l)[H,]  +  C  )t] 


-  exp(-k^[H2]t)} 


where  Ci  -  (lAykiCv) ) /k^  and  C2  -  The  term  C2  describes  the  re¬ 

moval  of  HF^  in  the  complete  absence  of  collislonal  quenching.  Within  the  ap¬ 
proximations  made  for  this  kinetic  scheme,^  the  rise  time  of  the  fluorescence 
Intensity  Tj.  is  directly  related  to  the  overall  rate  coefficient  k^  by  • 
(kjlH2])  and  the  decay  time  Tp  can  be  expressed  as  -  (k^(l)[H2]  +  02)*  A 
simplified  working  expression  for  the  fluorescence  intensity  is  given  by 


Eq.  (7): 


For  the  earlier  room  temperature  results,  this  analytic  equation  was  quantita¬ 
tively  tested  with  the  use  of  The  Aerospace  Corporation's  numerical  modeling 
code,  NEST^^.  In  that  study,  computer  modeling  calculations  indicated  that  as 


long  as  «  Tj,  then  is  a  good  approximation.  For  comparison 

with  the  data  of  the  present  study,  the  modeling  calculations  were  performed 
at  several  temperatures  between  200  and  1000  K.  Literature  values  were  used 
for  the  Einstein  coefficients,^^  the  relative  distributions  over  the  product 
vibrational  states, and  the  rate  coefficients  for  HF(v)  deactivation  by 
and  The  V-R,  T  deactivation  of  HF(v)  by  HF  was  assumed  to  occur 

by  single-quantum  processes  with  rate  coefficients  scaling  as  v  for  the 
higher  vibrational  levels. The  temperature  dependence  used  for  HF(1)-HF 
relaxation  was  taken  from  Ref.  5.  The  calculations  Included  losses  caused  by 
radiative  decay  but  not  by  diffusion. 

The  calculated  fluorescence  traces  were  fitted  to  [exp(-t/T^)  - 
exp(-t/T^)]  and  and  were  determined  In  the  same  manner  as  In  the 
analysis  of  the  actual  experimental  data.  Calculations  were  performed  for 
several  H2  concentrations  at  the  same  temperature,  and  the  values 
of  and  were  plotted  vs  H2  concentration.  It  was  found  that 

^  kj^[H2l  within  the  accuracy  of  the  curve  fitting  as  long  as  there  Is 
sufficient  H2  so  that  it  will  not  be  depleted  during  the  reaction.  The  rise- 
times  were  not  sensitive  to  changes  in  the  relative  product  distributions  or 
the  Einstein  coefficients. 

The  calculated  fluorescence  traces  decreased  from  their  maximum  to  about 
e”^  of  their  maximum  with  approximately  exponential  decay  rates.  A  plot  of 
decay  rates  at  700  K  vs  [H2]  had  a  slope  of  3.3  *  10^^  cm^/mol-sec  compared 
with  the  value  of  5.4  x  10^^  cm^/mol-sec  used  for  k^(l),  the  rate  coefficient 
for  HF(1)  deactivation  by  H2.  Thus,  the  H2  quenches  the  total  fluorescence  of 


Che  highly  excited  HF  with  a  quenching  rate  ~0.5S  tines  the  deactivation  rate 
of  HP(l)  by  H2. 

The  decay  times  obtained  fron  the  modeling  calculations  are  longer  than 
the  measured  values;  the  H2  quenching  rate  of  the  total  fluorescence  Is  about 
twice  as  fast  as  predicted  by  the  modeling  calculations.  An  experimental 
explanation  could  be  low-level  impurities  of  H2O  in  H2  although  care  has  been 
taken  Co  avoid  such  Impurities.  Ocher  effects,  such  as  diffusion  out  of  the 
observed  volume  or  deactivation  by  SF^  or  SF^  fragments,  would  not  be  propor¬ 
tional  to  the  H2  concentration  and  tfould  not  contribute  to  the  apparent 
quenching  rate  of  the  fluorescence  by  H2>  The  source  of  the  discrepancy  may 
also  be  the  computer  model.  The  modeling  calculations  were  performed  with 
single-quantum  processes  for  the  deactivation  of  HF(v)  by  HF  and  H2»  A  study 
by  Douglas  and  Moore^^  indicated  Chat  HF(4)  la  deactivated  by  HF  in  a  single¬ 
quantum  process.  Multi-quantum  deactivation  of  HF(v)  by  H2  would  bring  the 
modeling  calculations  Into  better  agreement  with  the  experimental  data  but 
would  not  completely  resolve  the  discrepancy. 

In  Fig.  3(a),  we  have  plotted  the  values  of  At^^/A[H2]  obtained  from  the 
experimental  data  and  from  the  computer  modeling  calculations  vs  1/T  (these 
latter  values  are  the  slopes  of  decay  rates  plotted  vs  [H2]  concentration 
at  a  constant  temperature).  For  comparison  with  these  data,  the  values  of 
^^(v  *  1)  are  also  plotted.  The  expected  value  of  At^^/A[H2]  is  k4(v  •  1)  if 
k^(v)  has  a  particular  v  dependence.^  In  Table  I  the  results  of  Che  modeling 
calculations  performed  at  five  different  temperatures  are  listed.  These 
results  Indicate  that  AT^^/kj^A[H2l  ■  0.99  ±  0.03,  which  Is  within  the 
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Fig.  3.  (a)  /AtH,]  vs  1/T;  Qpresent  experimental  data  obtained  from 

mixtures  of  Ar  and  H2;  0 numerical  modeling  results  (see  text); 

- values  of  k^(l)  used  In  the  numerical  modeling. 

(b)  At^  /AfDjl  vs  1/T; 0  present  experimental  data  obtained  from 
mixtures  of  Ar  and  D2;  0  numerical  modeling  results  (see  text); 
- values  of  k4(l)  used  In  the  numerical  modeling. 
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TABLE  I.  Nunerlcal  oodeling  results  for  representative 

conditions  of  5.4  x  nol/cnr  of  F  in  Ar  +  H2 

at  a  total  pressure  of  5.2  Torr. 


Temp  (K) 

AT“Vk^(l)A[H2l 

200 

1.00 

0.50 

295 

0.98 

0.53 

450 

0.98 

0.58 

700 

1.02 

0.60 

1000 

0.97 

0.55 

Average  0.99  ±  0.03 


Average  0.55  ±  0.05 


accuracy  of  the  graphical  data  reduction.  The  values  of  AT”^/k^(l)A{H2]  are 
close  to  0.55  ±  0.05  over  the  entire  temperature  range. 

Modeling  calculations  performed  for  the  F  +  D2  experimental  conditions 
indicated  that  »  k2(D2l  Is  a  good  approximation.  However,  the  discrepancy 
between  the  calculated  values  of  At^^/A[D21  and  the  experimental  values  shown 
In  Fig.  3(b)  Is  even  larger  than  the  discrepancy  for  the  F  +  H2  data  In 
Fig.  3(a).  The  F  +  D2  reaction  produces  DF(v)  v  <  4  compared  with  vibrational 
levels  V  <  3  produced  In  the  F  +  H2  reaction,  so  the  numerical  modeling  of  the 
^+02  experiments  Is  even  more  dependent  on  the  assumption  of  single-quantum 
transitions  and  the  v  dependence  of  the  DF  relaxation  rate  coefficients. 

The  multlphocon  dissociation  process  Is  affected  by  both  pressure  and 
temperature.  Our  experiments  are  performed  at  a  constant  Ar  buffer  pressure 
or  5  Torr.  Bado  and  van  den  Bergh^^  found  a  very  small  effect  of  Ar  pressure 
on  the  SF^  dissociation  yield  at  room  temperature.  We  are  not  concerned  with 
the  question  of  whether  the  multiphoton  dissociation  Is  colllslonless  In  our 
experiment.  The  low  Intensity  absorption  by  SF^  at  the  P(20)  CO2  laser  line 
has  been  studied  at  elevated  temperature  by  Nowak  and  Lyman. The  absorption 
coefficient  drops  rather  sharply  between  500  and  1000  K  as  a  result  of  red- 
shlftlng  of  the  spectrum.  The  multlphoton  dissociation  spectrum  Is  known  to 
be  red-shifted  from  the  absorption  spectrum,  primarily  for  this  reason.  Thus, 
the  optimum  P-branch  line  for  SFg  multlphoton  dissociation  certainly  shifts  to 
higher  J  as  temperature  Is  increased.  If  studies  at  T  >  800  K  are  pursued  In 
an  Improved  geometry  that  limits  background  Interference,  these  phenomena  may 
be  Important  experimental  considerations. 
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In  our  earlier  paper, ^  we  argued  that  the  F  atoms  were  translatlonally 
equilibrated  by  the  excess  of  Ar  to  H2(>S0  to  1)  before  reaction  with  H2«  In 
addition,  all  available  evidence  suggests  that  the  F  atoms  are  produced  with 
very  little  excess  translational  energy  from  the  multiple  photon  dissociation 
of  SFg.  The  question  of  electronic  equilibrium  was  not  considered.  In 
Sec.  V,  the  status  of  the  work  on  the  relative  reactivities  of 
F(^F3/2)  Is  reviewed.  Although  we  know  of  no  experimental  data  for  the  elec¬ 
tronic  quenching  of  exist  for  Cl(^Pj^2)  quenching. In 

the  case  of  Cl(2pj^2)  where  the  spin-orbit  splitting  Is  881  cm"^,  the  quench¬ 
ing  coefficient  by  Ar  Is  (1.1  ±  0.3)  x  10”^^  cm^ /molecule-sec.  Although  the 
quenching  of  f(  P1/2)  Ar  Is  expected  to  be  faster  since  the  spin-orbit 
splitting  Is  only  404  cm~^,  the  presence  of  these  two  thermally  accessible 
states  Is  cause  for  concern. 
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IV.  RESULTS 


With  the  use  of  the  Information  tabulated  In  Table  I,  we  can  ascertain 
whether  the  analytic  equations  (6)  and  (7)  are  valid  for  experimental 
fluorescence  traces.  The  slopes  of  the  experimental  vs  [H2]  and  [D2] 
curves  are  plotted  vs  1/T  In  Fig.  3  as  are  the  slopes  of  the  NEST- 
slmulated  x^^  vs  [H2]  and  [D2]  curves.  Removal  of  H2O  from  the  Ar  and  Ar  -f 
H2(D2)  flows  was  the  primary  cause  of  the  Improved  agreement  compared  with  our 
original  study.  Nevertheless,  there  Is  substantial  disagreement  between 
experiment  and  model.  No  experimental  data  are  available  for  HF(v  *  2  and  3) 
or  DF(v  ■  2,3,  and  A)  quenching  at  high  temperatures  by  H2  and  D2,  respec- 
tlvely.  Recent  measurements  at  T  ■  200  and  295  K  Indicate  that  the  vibra¬ 
tional  scaling  of  these  quenching  rates  Is  rather  Insensitive  to  tempera¬ 
ture.  As  a  first  approximation  In  the  modeling  calculations,  all  HF(v)  levels 
were  assigned  the  T-dependence  found  for  HF(1)  and  all  0F(v)  levels  that  found 
for  DF(1). 

Table  I  demonstrates  that  the  Inverse  rise  times  of  the  HF^  and  DF^ 
fluorescence  x^^  are  given  to  an  excellent  approximation  by  k^[H2l  and  k2[D2l, 
respectively,  at  all  temperatures  studied.  Plots  of  x^^  vs  [H2}  and  [D2]  are 
presented  In  Fig.  A.  The  resultant  values  of  kj^  and  k2  obtained  from  a  linear 
least  squares  fit  to  these  data  are  presented  In  Table  II  as  a  function  of 
temperature.  The  quoted  errors  In  the  individual  kj^  and  k2  measurements 
represent  the  sum  of  the  statistical  error  (lo)  of  the  least  squares  fits  In 
Fig.  A  and  a  systematic  error  of  ±15Z,  which  primarily  reflects  measurement 
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ID?)  X  109|mol/cm3| 


Inverse  rise  time  t  ^  as  a  function  of  teaperature  for 
(a)  HF'  as  a  function  of  [H2]  and  (b)  DF*  as  a  function 


TABLE  II •  Absolute  rate  coefficients  as  a  function  of  temperature 


F  +  H2  — HpT  +  H 

T(K) 

kj^(T),  cm^/mol-sec 

295 

(1.8  ±  0.3)  X  10^^ 

297 

(1.8  ±  0.4)  X  10^^ 

333 

(2.3  ±  0.6)  X  10^^ 

369 

(2.8  ±  0.5)  X  10^^ 

477 

(4.0  ±  1.0)  X  10^^ 

557 

(4.5  ±  1.0)  X  10^^ 

625 

(5.6  ±  1.6)  X  10^^ 

765 

(6.6  ±  1.3)  X  10^^ 

ko  X 

F  +  D2  -4.  DF'  +  D 

T(K) 

k2(T),  cm^/mol-sec 

295 

(9.5  ±  1.7)  X  lo’-^ 

301 

(8.6  ±  1.7)  X  10^^ 

332 

(1.25  ±  0.3)  X  10^^ 

374 

(1.0  ±  0.2)  X  10^^ 

401 

(1.3  ±  0.6)  X  10^^ 

473 

• 

00 

M- 

0 

• 

X 

0 

553 

(2.2  ±  0.5)  X  10^^ 

678 

(2.9  ±  0.7)  X  10^^ 
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errors  In  pressure,  flow  rates,  and  temperature.  These  data  for  and  k2  are 
plofed  vs  1/T  In  Fig.  5. 

Over  the  temperature  range  studied  (T  ■  295-795  K) ,  no  non-Arrhenius  be¬ 
havior  could  be  detected  within  the  experimental  uncertainty  (see  Sec.  V). 
Thus,  a  least  squares  fit  to  the  equation 

In  k  =  A  -  Eg/RT  (8) 

was  made  using  the  transformation  of  statistical  weights  described  by 

Cvetanovlc  et  al.^^»^^  Since  all  experimental  values  of  kj^  and  k2  have  a 

nearly  constant  fractional  error,  this  Is  equivalent  to  performing  an 

unweighted  linear  least  squares  of  In  k  versus  1/T.  These  Arrhenius  plots 

yield  the  following  values  for  k-j^  and  k2:  k^  =  (1.3  ±  0.25)  x  10^^ 

exp(-(1182  ±  100)/RT];  k2  »  (6.4  ±  2.2)  x  10^^  exp[-(1200  ±  142)/RT]; 

+  H  +  D  *  exp[(18  ±  250)/RT]. 

2  2 


Temperature  dependence  of  the  absolute  reaction  rates  (a) 
(this  work):  WH  (Ref.  6);  C  (Ref.  10);  MF  (Ref.  40); 
(b)  HBGM  (this  work);  WH  (Ref.  6);  C  (Ref.  10). 


V.  DISCUSSION 


Room  temperature  data  for  F  +  H2  and  F  4-  D2  ^.re  reviewed  In  Ref*  1. 
Additional  detail  is  given  in  the  more  comprehensive  reviews  of  Refs.  2 
through  5.  Tables  III  and  IV  contain  the  Arrhenius  parameters  obtained  for 
the  F  +  H2  and  F  +  D2  reactions,  respectively.  A  variety  of  experimental 
methods  has  been  employed;  Wurzberg  and  Houston^  also  used  the  infrared 
multiphoton  dissociation-infrared  fluorescence  technique  in  the  complementary 
temperature  regime  of  T  °  190-373  K.  In  this  section,  these  two  studies  are 
compared,  and  the  predicted  ratio  k^/k2  is  compared  to  very  precise 
measurements  made  by  Persky^^  and  by  Grant  and  Root.^^ 

Table  III  and  Fig.  5(a)  clearly  indicate  a  systematic  error  (~30Z) 
between  the  present  data  and  the  Wurzberg-Houston  studies  in  the  overlapping 
temperature  range.  Neither  study  has  revealed  the  cause  of  this 
discrepancy.  The  agreement  is  much  better  for  F  +  D2,  as  Indicated  in 
Table  IV  and  Fig.  5(b).  Both  studies  deviate  from  the  relative  rate  measure¬ 
ments  kj^/k2  given  by  Persky^^  and  by  Grant  and  Root;^^  however,  the  deviations 
are  rather  slight.  The  recent  low-temperature  study  by  Bulatov^^  gives  a 
ki/k2  value  in  excellent  agreement  with  Wurzberg  and  Houston.^ 

Figures  5(a)  and  5(b)  indicate  the  possibility  of  concave  curvature  to 
the  Arrhenius  plots  for  these  two  reactions.  The  exact  one-dimensional 

3S 

quantum  mechanical  calculations  by  Schatz,  Bowman,  and  Kuppermann  for  F  +  H2 
and  F  +  D2^^  predict  concave  non-Arrhenius  behavior  over  the  temperature  range 
200  to  1000  K  on  the  basis  of  substantially  different  temperature  dependencies 
for  reaction  into  the  various  product  vibrational  channels.  The  experimental 
data  of  Perry  and  Polanyl^^  over  the  same  temperature  range  do  not  agree  with 
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temperature-dependent  rate  coefficients 


this  prediction.  These  data  Indicate  extremely  minor  temperature  variation  of 
the  Important  vibratlonally-specif Ic  rate  ratios  k^(3)/kj(2)  for  F  +  H2»  and 
k2(4)/k2(3)  or  k2(2)/k2(3)  for  F  4-  02*  These  one-dlmenslonal  calculations 
appear  to  be  Inadequate  to  describe  the  details  of  the  temperature-dependent 
experimental  reaction  rates. 

Transition  state  theory  (TST)  calculations  performed  by  Cohen^®  Indicate 
that  curvature  of  this  magnitude  would  be  expected  for  either  bent  or  linear 
Intermediate  F-H-H  and  F-D-D  complexes,  and  sample  curves  have  been  super¬ 
imposed  In  Figs.  5(a)  and  5(b). 

These  TST  calculations  generally  do  not  deal  with  a  problem  that  has 
plagued  classical  trajectory  studies  of  this  system.  Truhlar  and  Muckerman 
and  Newton  noted  that  F(  P3/2)  reacts  with  H2  on  two  doubly-degenerate 
potential  surfaces  in  Cg  symmetry.  The  lower  (1  ^A')  correlates  with  ground 
state  products  [HF(^r)  +  H(^S))  while  the  upper  (^A")  correlates  with  excited- 
state  HF[HF(^n)  +  H(^S)]  and  will  be  nonreactive. 

Similarly  ^(^^1/2)  Interacts  with  H2  only  on  the  2  ^A'  surface  that  also 
correlates  with  electronically-excited  products.  In  the  absence  of  non- 
adlabatlc  transitions,  one-half  the  collisions  are  potentially  reactive  if  all 

7 

the  F  atoms  are  in  the  P3/2  state.  This  value  decreases  to  a  limiting  value 
of  1/3  as  T  Is  increased,  thus  populating  the  nonreactlve  ^Pi/2  l®vel.  At 
Intermediate  temperatures,  this  multiple  surface  coefficient  Is  described  by 

[exp(-A/kT)  +  2]"^  (9) 

where  A  Is  the  spin  orbit  splitting  In  the  F  atom  (404  cm”^).  If  F  +  H2  and 
F  +  D2  collisions  are  highly  adiabatic,  the  above  discussion  Indicates  a 
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contribution  to  a  convex  Arrhenius  plot.  Tully^^  calculated  contributions  to 
nonadlabatlc  behavior  arising  both  from  a  breakdown  in  the  Bom-Oppenhelmer 
approximation  and  from  spin-orbit  coupling.  The  final  result  Indicated  that 
F(  P1/2)  was  roughly  10%  as  reactive  with  H2  as  was  F(  Pj^2)*  Thus,  concern 
about  the  presence  of  the  imiltiple  potential  surfaces  seems  warranted,  and, 
based  on  Eq.  (9),  a  pre-exponential  factor  less  than  or  equal  to  one-third  to 
one-half  gas  kinetic  should  not  be  surprising.  In  general,  the  Wurzburg- 
Houston  study^  and  the  present  study  are  consistent  with  that  interpreta¬ 
tion.  Muckerman  and  Falst^^  recently  published  a  quasi-classical  trajectory 
study  for  F  +  H2  that  explicitly  accounts  for  the  multiple  surface 
coefficient.  The  predicted  pre-exponential  is  quite  close  to  that  found  by 
Wurzberg  and  Houston.^  The  potential  surface  employed  was  the  widely  accepted 
"Muckerman  5,"  which  has  a  barrier  height  of  ~1.1  kcal/mol;  the  resulting 
T-dependence  of  k^^  indicates  an  activation  energy  substantially  larger  than 
that  found  by  Wurzberg  and  Houston  or  the  present  study.  Calculations  on  a 
surface  with  a  slightly  smaller  barrier  would  provide  an  interesting  compari¬ 
son  between  theory  and  experiment. 

Over  the  temperature  range  studied  by  Wurzberg  and  Houston  and  ourselves, 
arguments  for  non-Arrhenius  behavior  are  ambiguous.  The  transition  state 
theory  calculations  by  Cohen^^  offer  a  plausible  fit  to  the  data;  however,  the 
transition  state  theory  methodology  is  sufficiently  flexible  to  allow  a  good 
fit  by  a  number  of  different  Initial  assumptions.  Certain  of  these  points 
have  been  recently  reviewed. The  importance-sampling  method  used  by 
Muckerman  and  Faist^^  should  allow  trajectory  calculations  of  sufficient  pre¬ 
cision  to  test  for  non-Arrhenius  behavior  in  the  200  to  1000  K  range.  We  do 
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not  believe  that  differences  In  the  apparent  of  ~250  cal/mol  should  be 
dismissed  as  experimental  artifacts  given  the  wide  temperature  range  studied. 

In  the  temperature  range  200  to  1000  K,  between  65  and  55Z  of  Che  energy 
of  the  F  -t-  H2  and  F  +  D2  reactions  (^32  kcal/mol)  appears  as  product 
vibration. Therefore,  35  to  45Z  of  this  energy  (~13  kcal/mol)  Is  available 
to  heat  Che  remaining  reactants  In  an  unbuffered  laser  mixture  where  the  heat 
capacity  C^(T)  >  6-7  cal/mol  K.  Adiabatic  temperature  rises  of  1000  to  2000  K 
are  thus  possible,  even  In  the  absence  of  vibrational  relaxation.  This 
problem  Is  magnified  in  the  H2'-F2-02  chain  laser  where  the  reaction 
H  +  F2“*’HF(v)  +  F  releases  a  substantial  fraction  of  its  exothermicity  Into 
translational  and  rotational  degrees  of  freedom. 


VI.  CONCLUSIONS 


The  Infrared  multiphoton  dissociation-infrared  chemiluminescence  method 
has  been  applied  to  the  measurement  of  absolute  rates  of  F  +  H2  and  F  -l-  D2  at 
elevated  temperatures.  The  validity  of  using  Infrared  chemiluminescence  from 
HF^  and  DF^  to  monitor  F-atom  disappearance  has  been  numerically  modeled  and 
verified.  The  results  are  In  substantial  agreement  (although  not  completely 
for  F  +  H2)  with  those  of  Wurzberg  and  Houston.^  Rather  than  force  a  single 
Arrhenius  equation  to  fit  both  the  Wurzberg-Houston  data  and  the  present  data, 
we  have  Indicated  that  small  differences  In  the  reported  activation  energies 
in  the  190  to  373  and  295  to  765  K  temperature  ranges  are  entirely  plausible. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operatlona  of  The  Aerospace  Corporation  is  conducting  exper* 
Imental  and  theoretical  investigations  necesssry  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  Bllltary  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-temperature  thermomechanlcs,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory;  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory;  Microelectronics,  GaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation^  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory;  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radiation 
froa  the  atmosphere,  density and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics.  Infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 


